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ABSTRACT 

Using data acquired as part of a unique Hubble Heritage imaging program of broadband colors of 
the interacting spiral system M51/NGC 5195, we have conducted a photometric study of the stellar 
associations across the entire disk of the galaxy in order to assess trends in size, luminosity, and local 
environment associated with recent star formation activity in the system. Starting with a sample of 
over 900 potential associations, we have produced color-magnitude and color-color diagrams for the 
120 associations that were deemed to be single-aged. It has been found that main sequence turnoffs 
are not evident for the vast majority of the stellar associations in our set, potentially due to the 
overlap of isochronal tracks at the high mass end of the main sequence, and the limited depth of our 
images at the distance of M51. In order to obtain ages for more of our sample, we produced model 
spectral energy distributions (SEDs) to fit to the data from the GALEXEV simple stellar population 
(SSP) models of Bruzual and Chariot (2003). These SEDs can be used to determine age, size, mass, 
metallicity, and dust content of each association via a simple \ 2 minimization to each association's 
B, V , and /-band fluxes. The derived association properties are mapped as a function of location, 
and recent trends in star formation history of the galaxy are explored in light of these results. This 
work is the first phase in a program that will compare these stellar systems with their environments 
using ultraviolet data from GALEX and infrared data from Spitzer, and ultimately we plan to apply 
the same stellar population mapping methodology to other nearby face-on spiral galaxies. 
Subject headings: catalog - Hertzsprung-Russell diagram - stars: early-type - galaxies: open clusters 
and associations: general - galaxies: photometry - galaxies: individual (M51) 



1. INTRODUCTION 

M51(also known as NGC 5194) is one of the nearest 
and largest grand design spiral galaxies. The galaxy is 
known to be undergoing an interaction with its compan- 
ion NGC 5195, and shows an active system of star forma- 
tion (SF), largely confined to two spiral arms. The mor- 
phology of NGC 5194 is listed by de Vaucouleurs 1976 
as a SA(s)bc(pec)III galaxy. The companion has been 
determined to be an SB0 (Sandage & Tammann 1987), 
and as an 10 (pec) by de Vaucouleurs, and its nucleus 
hosts a nonthermal source (Ford et al. 1985). Physical 
parameters adopted in this paper include a distance of 
8.4 ± 0.6 Mpc (Feldmeier, Ciardullo, & Jacoby 1997), a 
tilt angle of 20° and a positional angle of the major axis 
of 170° (Tully 1974). M51 is consistent with super-solar 
metallicity (Z ~ 2-3 Zq) and shows a weak metallic- 
ity gradient as a function of distance from the nucleus 
(Zaritsky et al. 1994). 

1.1. M51 Dynamics 

The M51 system is a very complicated and dynamic 
system. Most of the unique properties exhibited by the 
system appear, not surprisingly, to be directly related 
to the interaction the disk is experiencing with the com- 
panion. This interaction has been modelled extensively 
throughout the years, with increasing success in match- 
ing the present state of the M51/NGC 5195 system. 
Early models fit optical data well, but were unable to 
explain M51's extended HI tail (Toomre & Toomre 1972; 
Hernquist 1990; Rots et al. 1990). More recent attempts 
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fit HI maps well, but differ in the histories they paint for 
the M51 system. Salo & Laurikainen's (2000) N-body 
model suggests elliptical orbits and multiple passages of 
NGC 5195 through the disk of M51. In contrast, Wahde 
and Donner's (2001) model utilizing a genetic algorithm 
favors a hyperbolic and thus single-passage encounter, 
but fails to match the observed velocities in the HI tail. 
Theis and Spinneker (2003) combine N-body methods 
with a genetic algorithm, and achieve a good fit to both 
the HI intensity and velocity maps. Their results sug- 
gest a highly elliptic orbit with two recent disk passages 
at 400-500 Myr and 50-100 Myr ago, in agreement with 
Salo & Laurikainen. 

1.2. Star Clusters and Associations in M51 

In addition to its intrinsic properties, the uniqueness 
of the M51 system as a test bed for understanding the 
effects of galaxy interaction on star formation has led to 
many systematic studies of its stellar populations. Pre- 
vious HST studies of the compact star clusters in M51 
include Bik et al. 2003 and Bastian et al. 2005a,b, who 
investigated clusters near the center of M51, Lee et al. 
2005, who utilized a morphological cluster selection cri- 
teria across the whole of the galaxy, and Hwang & Lee 
2008, who detected over 2,000 clusters in the M51/NGC 
5195 system. 

Observations have shown that M51 exhibits a much 
bluer color than its companion overall, and the blue and 
yellow-red stellar populations within M51 are spatially 
separate structures (Zwicky 1955). The bluer star cluster 
colors exhibited in late- type galaxies such as M51 indi- 
cate recent star formation, as the O and B stars that 
produce this light are short-lived. In addition, inter- 
acting galaxies tend to contain very luminous and mas- 
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sive young star clusters, as compared to solitary galaxies 
(Whitmore 1999, and references therein). In this paper, 
and subsequent work, we intend to use the data from 
the ACS/WFC 10452 imaging survey of M51 to study in 
detail the structure, dynamics, stellar populations and 
ongoing star formation in this system. 

The high spatial resolution and broad wavelength cov- 
erage of the ACS /WFC provide unique data on the stel- 
lar populations in and around HII regions in M51. Our 
initial objective is to observe colors and magnitudes of 
stellar associations covering a range in luminosity and 
galactocentric radius within M51. These data can be 
used to build color-magnitude (CM) and color-color (CC) 
diagrams for the youngest stars, and thereby assess the 
fundamental parameters of the young stellar populations. 
At the distance of M51, we calculate that the limiting 
magnitude of the ACS- WFC camera (My =-3. 12) will 
restrict detection of CM turnoffs to only the youngest 
and brightest stellar populations, corresponding to an 
age of ~12 million years in the U-band. Thus, the stel- 
lar groupings which we will explore here are generally OB 
associations (~100 pc in size) or subgroups within OB 
associations (~50 pc) (Brown et al. 1999, Elmegreen et 
al. 2006). The stellar data can also be used to assess pos- 
sible systematic variations in the stellar content of these 
complexes as a function of position within the galaxy. 

The timescale between local compression and the peak 
in SF can be made based on the apparent separation of 
the arms delineated by the dust and HII regions respec- 
tively (Shu et al. 1972). This timescale is found to be a 
constant at 4.6 Myr for the southern arm (Elmegreen et 
al. 1992), and increasing as a function of radius for the 
northern arm, although the actual value for the northern 
arm timescale is less than that exhibited by the southern 
arm. By looking at the stellar components across the di- 
vide and their apparent ages we can learn how long the 
youngest stars reside where they were born, and thus 
track the position of the stars over time and when and 
where the SF event actually happened. In addition, we 
can use this stellar association sample to gauge the SF 
history as a function of position and determine the role 
the local environment plays in the rate and efficiency of 
the process. We already know that environment plays a 
key role in determining the stability of the local ISM to 
triggers or compressions that may induce gravitational 
collapse, and ultimately cause star formation. What is 
less clear is how that characteristic is set by the presence 
of a first generation of stars from an earlier SF event. 
Clearly, if a massive stellar cluster is formed, that clus- 
ter imparts not only strong UV radiation to the local 
ISM (as the resulting HII region is not typically radia- 
tion bounded), but also leaks strong continuum photons 
into the surrounding gas. The cluster can also contribute 
strong mechanical energy from its intense stellar winds 
that can sweep up and collect material into dense sheets 
or interfaces. These interfaces are ripe for the second 
generation of SF to ignite upon the incidence of a su- 
pernova blast wave (from the same short-lived massive 
stars), or from the nearby compression of material due 
to a variety of other possible triggers, such as cloud-cloud 
collisions, spiral density waves, etc.. 



2. OBSERVATIONS 

The datasets used in this study were obtained with 
the Hubble Advanced Camera for Surveys/ Wide Field 
Camera (ACS/WFC) as part of Hubble Space Telescope 
(HST) program 10452, which was completed in Jan- 
uary 2005 and became publicly available in April 2005 
(Mutchler et al. 2005). The datasets provide a six-tile 
image mosaic, with four dithered exposures taken for 
each of four filters. The 4 different imaging bands are 
B, V, and / broadbands (F435W, F555W and F814W 
respectively), and the narrowband redshiftcd Ha filter 
F658N. The standard archival data was retrieved "on- 
the-fiy" from the archive with standard ACS pipeline 
(CALACS) processing, including bias, dark, and flat- 
field corrections. The Ha exposures are of sufficient 
depth to reach surface brightnesses of at least 10~ 16 
ergs/cm 2 /s/arcsec 2 to allow us to probe the faintest HII 
regions. 

An ACS/WFC pixel subtends 2.0 pc at the distance 
of M51 and its companion NGC 5195. This means that 
a single pixel in the image could contain one or multiple 
stars, and thus there is the possibility of more than one 
star being sampled by each pixel. If two or more stars fall 
on the same pixel, this will still yield a non-symmetric 
wing in the PSF that may be detected in the surrounding 
pixels. Also, the total luminosity of a pixel is a good 
determinant of the number of stars in a single pixel or 
source. If the intrinsic brightness of a pixel or source 
is greater than the maximum luminosity expected for a 
single star then it can be assumed that the pixel contains 
multiple stars. On our CM diagrams, luminosities at 
or above My— ~-9 are assumed to be multiple stars or 
compact clusters. 

2.1. Pointings 

Using the Wide Field Channel (WFC) of the Advanced 
Camera for Surveys (ACS), a 2x3 mosaic (6 pointings, 
4 orbits per pointing) was obtained in 4 filters (Johnson 
B,V,I, and Ha) at a telescope orient of ~270°. Details 
of these exposures along with the central coordinates of 
the 6 pointings or "tiles" in the 2x3 mosaic are available 
in Mutchler et al. 2005. Each pointing provides 4 orbits 
of integration providing the limiting magnitudes listed in 
Table 1. The 4 exposures within each tile were dithered: 
a small sub-pixel dither (2.5x1.5 pixels), and a larger 
dither which spans the interchip gap (5x60 pixels). 

3. DATA REDUCTION 

3.1. Processing 

The initial raw data was processed using the stan- 
dard Hubble pipeline, through the IRAF task CALACS 
(Mutchler et al. 2005). CALACS performs the correc- 
tions for bias, dark, and flat fielding. The resulting im- 
ages were then cosmic-ray rejected by combining two im- 
ages of the same exposure time and gain, discarding the 
saturated cosmic ray pixels from one image, and replac- 
ing the "bad" pixel with the good signal from the second 
image. Distortion effects from the camera's off-axis posi- 
tion within the observatory are smaller than the photom- 
etry aperture used, and thus are inconsequential for this 
application (but would be relevant if exact astrometry of 
individual stars was needed). 
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3.2. Photometry and Initial Stellar Association 
Selection 

The associations included in this catalog were chosen 
to be relatively isolated from neighboring groups of stars 
in order to produce a single-aged CM diagram. After 
processing with CALACS and cosmic-ray rejection, the 
V-band image was inspected by eye for isolated stellar 
groups from which CM and CC diagrams could be de- 
rived. Initial sample selection was performed by visually 
scanning each image for groupings of stars/sources that 
could be associated with one another. A total of 969 
stellar associations were selected in this way for subse- 
quent analysis. For each association, a CM diagram was 
produced to assist in determining if the group of stars 
in question were indeed coeval. A more detailed descrip- 
tion of this method of analysis is described in follow- 
ing sections. Photometry was performed on these group- 
ings of stars using the GSFC IDL ASTROLIB package 
IDLPHOT. This package selects the individual stars and 
outputs the flux from each star, as measured within an 
aperture with a two pixel radius, in the B, V, and I 
broadband filters and the Ha narrowband filter. Two 
pixels is the approximate full width at half maximum 
(FWHM) of point sources in these images, and this aper- 
ture size was chosen to help exclude unresolved compact 
clusters from our source selection, although selection of 
the most compact clusters along with the stars is un- 
avoidable. Thus the "sources" plotted include stars, un- 
resolved multiple star systems, and some compact clus- 
ters. Background maps for each tile, chip and imaging 
band were produced by iteratively replacing any pixel 
greater than 1.3 times the median of the surrounding 
12x12 pixel box with the median value, until all stellar 
sources were masked out by the median value. These 
background maps were then subtracted from the images 
to properly correct for the local sky level. The fluxes 
derived from our initial methods were then converted 
to absolute Vega magnitudes by dividing by the expo- 
sure time, multiplying by the spectral flux density per 
unit wavelength that would generate 1 count / sec (header 
value PHOTFLAM), and then translated to an ST mag- 
nitude via the equation m = -2.5 logi (F) + PHOTZPT 
(where PHOTZPT=-21.10 for the ST magnitude scale). 
These magnitudes were converted to Vega magnitudes 
and compared to magnitudes calculated via the method 
detailed in Sirianni et al. (2005). The magnitudes cal- 
culated in these two methods were identical to within a 
thousandth of a magnitude. Finally, a CM diagram of 
absolute magnitude in the V-band versus (B — V) and 
a CC diagram of (B — V) versus (V — I) were produced 
for each stellar association. 

3.3. Extinction/Reddening and Redshift Corrections 

Photons collected from the M51/NGC 5195 system ex- 
perience reddening and extinction due to dust internal to 
M51 and Milky Way dust, as well as redshift from the 
radial component of the galaxy's proper motion. All of 
these factors must be carefully corrected to insure ac- 
curate photometry. The redshift of M51 is z—0. 00154, 
which corresponds to shifts of 6.62A, 8.23A, 12. 8A, and 
10.1 A, for the central wavelengths of the B, V, I, and 
Ha filters, respectively (Scoville et al. 2001). These cor- 
rections were applied to all wavelengths. Extinction due 



to dust poses a more complex problem, as decomposing 
internal and Galactic dust is extremely difficult without 
data across the whole spectral range. Also, most authors 
utilize the Galactic extinction corrections from Table 6 of 
Schlcgel et al. 1998, but these values are not ideal for our 
purposes since they were acquired by convolving the SED 
of an elliptical galaxy with the specified filter through- 
put curves. Assuming the SED of an elliptical galaxy 
neglects the blue continuum from young stars present in 
a spiral galaxy, and would introduce an ~5% error in 
our reddening corrections (Schlegel 2008, personal cor- 
respondence). For the present analysis, we treat the in- 
ternal and Galactic extinction separately, and follow the 
precedent of other authors (Lee et al. 2005, Mora et al. 
2009, Ubeda et al. 2007) and use the Schlegel et al. 1998 
values to correct for extinction due to dust in the Milky 
Way, in spite of the 5% error. Future work is planned to 
explore the possibility of avoiding this error by produc- 
ing de-reddened SSP SEDs from the GALAXEV models 
and fitting the reddening values to the data available 
from Hubble, GALEX, and Spitzer for each stellar as- 
sociation. This method would bypass the necessity of 
decomposing internal and Galactic reddening, and avoid 
the error imposed by assuming the SED of an elliptical 
galaxy. For the present study, we correct for Galactic ex- 
tinction using equations 1, 3a, and 3b in Cardelli, Clay- 
ton, and Mathis 1989, with i?y=3.1 and Ay MW =0.117 
for the direction of M51 (Cardelli et al. 1989, Schlegel et 
al. 1998). 

The intrinsically dusty environment of M51 is known 
to hamper studies of both the ionized gas and the HII 
regions themselves. The overall apparent structure and 
excitation state of gas around massive ionizing stars is 
strongly affected by the presence of the dust and may 
modify the continuum of photons that escape to ion- 
ize and energize the surrounding interstellar medium 
(ISM). The effect these modifications have on the dynam- 
ics and susceptibility of the ISM to subsequent second- 
generation triggering or propagation of star formation 
make this galaxy an interesting candidate for the study 
of star formation as a global mechanism for distributing 
and regulating the formation of the next generations of 
stars. To account for the effects of dust internal to M51, 
we assume a global extinction correction of Ay M51 =0.3 
mags, corresponding to an reddening correction of E(B- 
V)=0.097, using Ry=3.1 (Cardelli et al. 1989). This 
value was found to be an average of the best fit for the 
majority of the associations studied, if Ay M51 is allowed 
to vary between 0-1 mags, in increments of 0.1 mags. 
This value is in good agreement with with Hwang & Lee 
2008, who use E(B-V)=0.1 (A Vmbi =0.31 mags), and is 
less than the maximum value in the bulge, which is E(B- 
V)=0.25 (Avm 51 =0.775) (Lamers et al. 2002). Ideally, it 
would be best to leave the extinction correction as a free 
parameter in the SED fitting, since dust content varies 
noticeably over the disk of M51, and thus a single correc- 
tion value for the entire galaxy can only be an approx- 
imate correction. Unfortunately, with only BVI fluxes, 
we are unable to discern fluctuations in the local dust 
environment and stellar association ages simultaneously 
Visual inspection of the CM and CC diagrams of the 
120 associations in our set show that Ay M51 lies between 
and 1 magnitudes for most of the associations. We 
choose an intermediate value of 0.3 magnitudes, know- 
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ing that some associations will be slightly over- or under- 
corrected. Once redshift corrections and internal and 
Galactic reddening values. 

3.4. Stellar Association Selection and Age 
Determination 

The use of a CM diagram of My versus (B-V) was origi- 
nally chosen to allow direct comparison with the Sandage 
age estimation model (Sandage 1957), in order to test 
the viability of determining stellar association ages via 
their main sequence (MS) turnoff magnitude, given the 
distance of M51 and depth of HST ACS/WFC images. 
These colors were chosen over other available bands to 
maximize the depth of the sample, as B and V have 
fainter limiting magnitudes than I or Ha. The Sandage 
model is based on MS turnoff ages of 10 galactic clus- 
ters and one globular cluster. These clusters provide an 
age-MS turnoff magnitude relation, which we derived via 
a simple empirical interpolation of Sandage's Figure 1. 
For each association in our sample, a selection box of the 
appropriate size was chosen to produce the narrowest 
sequence and thus the tightest fit to the Sandage model 
when it was overlaid on the resulting CM diagram. Dense 
stellar regions which did not closely fit the model were 
not included in the list of potential associations, as they 
are most likely blends of multiple cluster. These groups 
may be spatially related, but are not necessarily coeval, 
making determination of a single age impossible. 

More modern treatments of this approach involve the 
use of model isochrones to accurately predict the appar- 
ent locus of stellar points on such a CM diagram for 
known differences in metallicity, age, and cluster size. 
However the underlying treatment is still sound - to use 
the distribution of member stars in an association to de- 
termine an estimate for the apparent age of the grouping 
based on the turnoff exhibited. Complications arise in 
the application of the method due to blends of multi- 
ple stars and/or compact clusters being detected as a 
single source. Plotting sources in both color-magnitude 
and color-color space helps to separate stars from com- 
pact clusters, and to weed out groups of stars that are 
not massive enough to be accurately fit with the Sim- 
ple Stellar Populations models of BC03. The M51/NGC 
5195 system is at the very limit of the HST/ACS cam- 
era's ability to resolve individual stars, and the corre- 
sponding brightness limit due to the exposure times used 
narrows the ages we could detect via this method to only 
the most recent of SF events. Consequently, this study 
is only sensitive to the bluest of associations containing 
the most massive, and therefore the shortest-lived, stars. 

Even if blending or resolution limits were not an is- 
sue, the fact that we are working at the blue end of the 
MS also constrains the resolution in magnitude space 
that can be used to resolve turnoffs on CM diagrams. 
At the top end of the MS, the sequence is almost ver- 
tical, meaning that the locus of turnoff points is also 
almost vertical. The degree to which we can measure 
small changes in brightness therefore limits the degree 
to which we can pull out departures from the MS due 
to evolutionary turnoff, as opposed to a simple intrin- 
sic spread in brightness due to population and possibly 
internal reddening etc. There is an intrinsic scatter in ap- 
parent brightness for any collection of stars in a cluster 
or association, and it is the art of distinguishing intrinsic 



scatter from changes due to stellar evolution that is at 
the heart of this approach, and is also the most difficult 
part of the technique. The limiting magnitude of the 
ACS observations translates into an inability to detect 
main sequence turnoff ages any older than ~12 million 
years. This is a flaw in the approach, but given the dis- 
tance to the galaxy, this is at the boundary of what we 
can do observationally at this time (Rey et al. 2007). 

To help make the final selection of single- 
aged associations to be included in the sam- 
ple^ Padova stellar isochrones (Girardi 2006, 
http:/ /pleiadi. pd.astro.it/isoc_photsys.02/isoc_acs_wfc/) 
and a BC03 stellar population evolutionary track were 
overlaid on the CM and CC diagrams. Figure 1 shows 
these plots for three example associations from our set 
that span the range of ages found (4 to 610 Myr). One 
must carefully assess the distributions of sources for 
each association on these two plots to be certian that the 
stars/compact clusters in it are in fact coeval. Groups 
of stars that did not lie along a single isochrone (within 
the scatter of the data) , whose total magnitude did not 
align with the BC03 model SED, or whose CM and CC 
diagrams indicated discrepant ages were excluded from 
the set. This second iteration of the selection process 
yielded a final list of 120 associations. Figure 1 shows 
the CM and CC diagrams for three representative stellar 
associations from this set, alongside the image subraster 
that generated it. 

Since M51 is at the limits of our ability to resolve indi- 
vidual stars and measure MS turnoff ages, we decided to 
employ a modern refinement - the use of stellar popula- 
tion models in tandem with SED fitting to allow the pre- 
diction and subsequent measurement of stellar associa- 
tion age. In order to fit model SEDs to the selected stellar 
associations, we have written a series of IDL scripts that 
calculate the flux from each of the associations, produce 
model stellar population SEDs from the GALAXEV li- 
brary, calibrate these models for comparison to the data, 
and perform a x 2 fit between the observations and the 
model SEDs. The SED with the minimum \ 2 value when 
compared to the data is the best-fit model. 

The GALAXEV library allows for tuning of the in- 
put parameters in order to best match the stellar pop- 
ulation of interest, in our case stellar associations in 
the M51/NGC 5195 system. For associations in M51 
we choose the Padova 1994 stellar evolutionary tracks 
(BC03, Alongi et al. 1993, Bressan et al. 1993, Fagotto 
et al. 1994a, b, Girardi et al. 1996), a Salpeter initial 
mass function of (;(log w (M)) cx M -1 - 35 (Salpeter 1955), 
and an instantaneous burst SFR of ip(t) = 1M® ■ S(t) 
(BC03), which best replicates the scenario of triggered 
star formation. We allow for a range of models between 
the ages of 0.1 Myr to 1 Gyr (this range should include 
the two most recent perigalacticon passages of the two 
galaxies) (Salo & Laurikainen 2000), extinction due to 
dust in M51 of __y=0.3 magnitudes, and a metallicity 
of 2.5__0 (Zaritsky et al. 1994). The mass of each asso- 
ciation is quantitatively calculated during the \ 2 fitting 
process, details of which are described below. 

To perform model SED fitting, a library of BC03 
Simple Stellar Population models was produced, and 
the model stellar association SEDs were created by 
converting the GALAXEV output into fluxes that 
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Fig. 1. — CM and CC diagrams, and V-band images for three example stellar associations from our set. The blue-purple dotted lines are 
Padova stellar isochroncs for stellar masses from 0.15 to 66.56 MSun, metallicity Z=0.019, and ages are 1 Myr (lightest blue), 6 Myr, 10 Myr, 
20 Myr, and 50 Myr (purple) (Girardi 2006, http://pleiadi.pd.astro.it/isoc_photsys.02/isoc_acs_wfc/). The red line in the left hand panel 
is a BC03 model evolutionary track for an instantaneous burst stellar population ot metallicity £=0.02. Datapoints have been corrected 
for Milky Way dust, and A^ M51 =0.3 mags has been added to all models plotted. The two large dots represent the total association flux 
from summing the sources (open circle), and the total flux from within a large circular aperture centered on the association(closed circle). 
Ages for these three associations are, from top to bottom, 4 Myr, 160 Myr, and 610 Myr. 
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can be compared directly to fluxes calculated from 
the Hubble data. The 'library of BC03 SSP mod- 
els are created by running two scripts available in the 
GALAXEV package: CSP_GALAXEV and GALAX- 
EVPL. CSP_GALAXEV evolves the composite stellar 
populations spectrally for a constant metallicity using 
equation (1) of Bruzual & Chariot 2003, and outputs a 
binary file of the evolved spectral evolutionary model. 
The program GALAXEVPL then reads in the binary 
file from CSP_GALAXEV, and outputs ASCII hies of 
SEDs at the ages requested by the user. This ASCII file 
contains wavelength (column 1) and specific luminosity 
density (£\ in the Papovich 2001), or luminosity per unit 
wavelength, in units of L Q ■ A -1 (column 2). In order 
to convert from this number to flux (F v ), we followed 
the procedure outlined in Section 4.2 of Papovich et al. 
2001, using 

"~ cM tot W 
to calculate the luminosity per unit frequency, and 



F v = (1 + z)L v 



And 2 



(2) 



to calculate the flux per unit frequency, where Ao is the 
redshifted wavelength, l\ is the specific luminosity den- 
sity, M tot is the total mass of the stellar population at a 
particular timcstcp, A\^mw is the extinction due to dust 
in the Milky Way (Cardelli et al. 1989, Schlegel et al. 
1998), A\.M5i is the extinction due to dust internal to 
M51, d is the distance to the galaxy, and all units are 
in the CGS system. F v calculated from equations (1) 
and (2) is for a 1 M Q stellar population, and has units 
of erg • cm~ 2 • s -1 • Hz -1 . The quantity M to t is ob- 
tained by extracting the value of Mgalaxy in column 9 of 
the *.4color ASCII file output from GALAXEVPL which 
corresponds to the desired age to be modelled, and is the 
same quantity as m*(t) in Papovich et al. 2001. 

We then convolved each model SED with the B, V, I, 
and Ha filter throughput curves to get the flux per unit 
frequency in each filter, /„, using 



I 



-dv 



f^du 



(3) 



where r„ is the total system throughput curve, including 
all transmission and detection effects of the telescope, 
camera, filter, and CCD chip (Fukugita et al. 1996). 

The mass of each stellar association is calculated using 
the fact that the BC03 model SEDs are normalized to a 
1 Mq stellar population. The mass which produces the 
smallest x 2 value between the data and each model SED 
is calculated by minimizing the following equation with 
respect to M and then solving for M, then plugging this 
mass back into the equation to calculate the minimum 
X 2 for each model. 



X 



E 



(4) 



where the summation is over i filters, f° b f is the observed 
flux in filter i, M is the best-fit mass for a particular stel- 
lar association, f v ^ is the model SED flux in the filter i, 



and Oi is the error in the observed flux for filter i. The 
best-fit model stellar association SED was then chosen 
from the array of model SEDS by finding the minimum 
X 2 value between model and data in the B, V, and / 
bands. It should be noted that it is not valid to use Ha 
flux directly as a fitting parameter for a stellar popula- 
tion, since Ha is emitted and absorbed both from the 
stellar photosphere and the surrounding gas. 

3.5. Other Parameters (Stellar Association Size, 
Position) 

The size of each association was estimated by deter- 
mining the spatial extent of the full width at half maxi- 
mum (FWHM) of the flux density profile of stars. Specif- 
ically, stars/sources were counted in radial bins from the 
association center, the flux from sources in each bin were 
totaled, and the standard deviation of the spatial flux 
density distribution was computed. The full width at 
half maximum was then computed assuming that the dis- 
tribution was Gaussian. 

The position of each stellar association in equato- 
rial sky coordinates was calculated using the WCSTools 
package (Mink 1997, 1999), which allows for the non- 
linear conversion from pixel coordinates in the undriz- 
zled images output from CALACS to right ascension and 
declination in the appropriate epoch. Associations were 
named according to the tile of the mosaic they reside in, 
and the lower left and upper right corners of the bound- 
ing box for that association, in pixels. 

4. RESULTS 

4.1. Properties of Isolated Stellar Associations in M51 

Results presented here are from SED fitting of the ages, 
masses, metallicities, and internal extinctions for each 
stellar association. Ages were also assessed by attempt- 
ing to pin down the location of the main sequence turnoff 
point on the color-magnitude diagram, but we found that 
turnoffs were not apparent for the vast majority of asso- 
ciations (see Discussion). Table 1 lists the derived prop- 
erties for the 120 associations selected from our sample 
using the criteria and methods laid out above. The iden- 
tifier for each association is based on a 4-digit number 
indicating which image in the mosaic the association is 
located in, and the x,y pixel coordinates of the lower left 
and upper right corners of the selection box used for se- 
lection of that association. The location of the central 
pixel in Epoch 2000.0 coordinates and BVI magnitudes 
in the Vega system are listed for each association. 

It should be emphasized that this is the subset of as- 
sociations from the full sample whose CM diagrams in- 
dicated that they were likely to be single-aged, allowing 
indisputable measurement of the properties listed, within 
the limitations of the method. 

Figure 2 shows a map of the galaxy with various anno- 
tations derived from the measurement exercise described 
above, and a number density histogram of the stellar as- 
sociation size sizes. On this diagram each of the associa- 
tions has been pictorially represented by a colored circle. 
The size of the circle reflects the relative size of the asso- 
ciation in physical extent. The color of the circle depicts 
the best-fit age of the association in question - with red 
associations being the oldest at ^600 million years, rang- 
ing down through the rainbow to the youngest ages found 
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at 4 million years, in blue. Figure 3 shows a similar map 
of the galaxy, where the size of the circle corresponds 
to the best-fit mass calculated quantitatively from the 
SED model, instead of the association size. Masses found 
range from -1.1 x 10 3 -7.4x 10 6 M . The contents of Ta- 
ble 1 are therefore pictorially represented in these two 
galaxy maps. 

5. DISCUSSION 

5.1. Age, Size, Mass, and Dust Distribution with 
Position 

Examination of Figure 2 (top) shows that the major- 
ity of our association sample is resident on one or the 
other of the two primary spiral arms in the galaxy, or 
on spurs coming off the spiral arms, with few found in 
the interarm region. This is not a surprising result, but 
considering that associations in our set are defined to 
be single-aged groupings, this makes a statement about 
the dominant outcome of SF in these regions due to the 
recent SF episodes active in the disk. The lack of single- 
aged groupings in the inter-arm region may indicate that 
stellar associations disperse relatively quickly after for- 
mation. There is also a noticeable concentration of asso- 
ciations along the northern arm, the arm closest to NGC 
5195, indicating an enhancement in the most recent SF 
activity in that area. The distribution of stellar associa- 
tion ages along this arm, coupled with the relative spatial 
positions of M51 and NGC 5195 indicates that gas is ac- 
tively streaming along the arm and forming stars in its 
wake. Clearly, our association selection criteria limits 
the sample to the most dominant SF regions, since we 
require the associations to be well defined and well iso- 
lated from neighboring associations or dense collections 
of field stars. Figure 2 shows no noticeable trend between 
global position and association age, but does show some 
local trends indicating the possibility of secondary and 
tertiary SF events. There is also no apparent correlation 
between association age and distance from the galactic 
nucleus measured along each spiral arm. The sequenc- 
ing of stellar associations spatially with age, which one 
would expect from SF propagating from a trigger, are 
apparent in some instances, particularly on the bend of 
the Northern spiral arm closest to the companion, but 
is not widespread across the whole of the galaxy. The 
cycling time nesessary to go from one generation of stars 
to the next is —5 Myr, meaning that many cycles of 
triggered star formation could have occured in this time 
range (4 to 610 Myr). It could be the case that while we 
see the sequential propagation of triggered star formation 
on star cluster scales (1-5 pc), this propagation may not 
occur on the larger stellar association and star complex 
scales (—50-200 pc) investigated in this study. Another 
interesting result is that we see no apparent surge in the 
local SFR at either of the two "breaks" in the spiral 
arms, where the gravitational effects of the interaction 
has warped the spiral structure. One might expect an in- 
creased SFR in these regions, but the lack of associations 
found here indicates that conditions are not conducive to 
the condensation of gas into stars. It is likely that the 
dynamical shear at these two "breaks" disrupts the col- 
lapse of H2 and HI clouds into stars. The size histogram 
shows a peak at ~25 pc, which may be due to the mixing 
or collapse length in the instellar medium dictating the 



dominant size of the fragment of the molecular gas cloud 
that will collapse into a stellar association. 

In Figure 3, we see that there is a tendency towards 
more massive stellar associations in the northeastern side 
of the galaxy, particularly near the nucleus and on the 
northern spiral arm. This trend hints at the fact that 
stellar association mass may be more sensitive to triggers 
such as galaxy interactions and spiral density waves than 
age or size of the association. 

The main sequence turnoff age determination method 
proved to have several major limitations. The vertical 
and horizontal location of any apparent MS turnoff is 
limited by the vertical slope of the MS at this location in 
the CM diagram, and the horizontal scatter of member 
stars in the association. Since our limiting magnitude of 
My=-3.12 defines that we are only seeing the very upper 
end of the MS, the OB stellar population. Typical life- 
times of the members of this population are really quite 
short, and indeed are much shorter than the 12 million 
year age cutoff that the M„=-3.12 should represent. As 
such many members of the associations we are study- 
ing never achieved turnoff, but instead probably went 
supernova. Older associations would, as a direct exten- 
sion of this principle, be composed of much fainter stars 
and would therefore fall below both our limiting mag- 
nitude and selection criteria for analysis as part of our 
sample, and not be resolved into individual populations. 
Analysis of turnoff tracks based on stellar evolution the- 
ory also indicates that those B stars that do survive this 
long follow much smaller turnoff tracks than we expect 
to see from lower mass stars further down the MS. All 
of these factors combine to diminish the size of the effect 
we are trying to measure and therefore act to reduce the 
final sample of associations from which we could draw 
conclusions that could withstand scrutiny. Technically, 
the main sequence turnoff method should allow detection 
of association MS turnoffs up to 12 million years of age, 
but no associations with a turnoff at this brightness were 
found, because of the confusion limit encountered when 
trying to verify the integrity of a turnoff right at the 
vertical magnitude limit of the observations on diagrams 
like the one depicted in Figure 1. 

5.2. The Recent Star Formation History of M51 

Using all the results presented in this paper, we can 
make several rather broad statements about the recent 
SF history of M51. It is clear from Figure 2 that the 
interaction with NGC 5195 has generated many large, 
well-isolated young stellar associations across the north- 
ern arm of M51, many more than are apparent along 
the southern arm. This implies an enhancement in the 
SF on that side of the galaxy, but does not represent a 
system- wide elevation in the SF rate (SFR). Within the 
limitations of the measurements we have been able to 
make as part of this study, we cannot see any systematic 
age sequencing of SF as a function of location. 

This is somewhat disappointing since one would expect 
to see evidence of age sequencing with location as a di- 
rect result of SF propagation in the disk. Why do we not 
see such evidence? Aside from the limits of the data, it 
is very possible that the ISM in the galaxy is sufficiently 
well-mixed and in an increased dynamic state that it is 
possible to trigger SF in any given location far easier than 
it might be with a more quiescent ISM. We are encour- 
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Fig. 2. — Top: A map of M51 indicating the age, size and location of each of the associations in our 120-association sample. For each 
association the relative size is represented by the size of the circle used, the age is represented by the color of the circle, and the location 
is presented as an overlay on the full field ACS V-band image itself. Bottom: A histogram of number of associations per size bin. 
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Fig. 3. — A map of M51 indicating the age, mass and location of the 120 associations in our sample. For each association the mass is 
represented by the size of the circle used, the age is represented by the color of the circle, and the location is presented as an overlay on 
the ACS V-band image. 
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aged that our sample reflects our original expectations 
concerning enhanced SF on the northern arm, it would 
however appear that that SF locations are much more 
random in nature. 

Indeed, the quiescent SF state reflected by the south- 
ern arm of the galaxy seems to represent a very low level 
of active SF over the last 5 million years, when compared 
to the northern arm. As noted in the Introduction, the 
expected timescale between the passage of the compres- 
sion wave in the spiral arm and the formation of stars is 
just less than 5 million years, and so this would suggest 
that these results are commensurate with this view, but 
that the size of the compression apparent in the northern 
arm, and therefore the SFR this represents, is dramat- 
ically elevated, and would result in a different Schmidt 
law. These results explicitly depict a large variation in 
the SFR active in the disk of M51 from one side to the 
other apparently as a direct result of the passage of NGC 
5195. 

5.3. Next Steps in the Analysis of the System 

While the human eye is an ideal detector for picking 
up patterns, the human brain is prone to bias. Thus our 
method of visual inspection for association selection and 
turnoff magnitude is prone to the bias of the observer. 
In order to reduce this effect, we plan to supplement 
our visual inspection methods by using an automated 
selection algorithm, to produce an initial set of candidate 
associations and ages, and then using visual inspection 
of each candidate to eliminate false positives and assure 
accurate age fits. 

Also, it has been shown that fitting model SEDs to 
optical broadband colors alone can be inaccurate, in that 
it is difficult to differentiate between the reddening effects 
of age, dust, and metallicity (Cervino & Luridiana 2006, 
Anders et al. 2004, Piskunov 2009). As such, in the next 
paper we will incoorporatc UV data from GALEX and 
IR data from Spitzer, in order to better differentiate the 
reddening effects of stellar population age, dust content, 
and metallicity. 

In the future, we plan to apply these same techniques 
to stellar associations in M83 and NGC 4214, using 
the multiwavelength images of this galaxy taken with 
the new Hubble Space Telescope Wide Field Camera 3 
(WFC3) . The wide range of wavelength sensitivity avail- 
able with WFC3 will allow for better differentiation be- 
tween the reddening effects of age, dust, and metallic- 
ity, and thus more accurate dating of stellar associations. 
Findings from the present study of M51 will be compared 
and contrasted with those from the non-interacting spiral 
galaxy M83. 

6. CONCLUSIONS 

We have taken the first high-reslution, multiwave- 
length, system-wide imaging dataset of the interacting 
spiral system M51/NGC 5195 and conducted a photo- 
metric study of the stellar associations in the system in 
an effort to measure and characterize the SF in the sys- 
tem, using the associations as age tracers. We identified 
969 stellar groupings of interest, and from this set drew 
a smaller sample of 120 associations that were deemed 
potential candidates for further analysis. Those associ- 
ations which were isolated enough to be unequivocally 
regarded as coeval and largely unpolluted by field stars 



were included in a final sample explored in this paper. 
This subset was extracted and presented as the best rep- 
resentees the properties of the stellar association popu- 
lation resident in M51. 

To summarize, from the selection, photmetry and 
SED fitting of 120 single- aged stellar associations in the 
M51/NGC 5195 system, we conclude the following main 
points: 

1. Visual inspection of the CM and CC diagrams for 
each potential stellar association is an effective way 
to weed out groups of star that are not coeval. We 
find that the vast majority of the single-aged stellar 
associations selected in this manner lie along one 
of the spiral arms on M51, or on spurs coming off 
the spiral arms. 

2. Very few single-aged associations are found in the 
interarm regions, suggesting that either the condi- 
tions in these inter-arm regions are not ripe for the 
formation of stellar associations, or that stellar as- 
sociations fromed here disperse quickly after their 
birth. 

3. There is an enhancement in the number of single- 
aged stellar associations in the northern arm spiral 
arm closest to the companion, most likely triggered 
by the intereaction of M51 and NGC 5195. 

4. No correlation between position within the galaxy 
and age of the association was found. There is also 
not an apparent relationship between the physical 
size of the association and position. 

5. We see few instances of sequential formation of stel- 
lar associations, which suggests that this type of 
triggering may only occur on scales smaller than 
investigated in this study (i.e. smaller ~5 pc, on 
the scale of compact clusters). 

6. Single-aged associations in this system range in size 
from 5 to 90 parsecs, with the peak in the size- 
number density histogram occurring at an average 
radius of ^25 parsecs, hinting that this may be the 
dominant size of the scale length for collapse within 
the molecular clouds that stellar associations form 
from. 

7. There is a slight trend for more massive associ- 
ations to to exist perferentially in the region of 
M51 closest to NGC 5195. This indicates that the 
galaxy interaction facilitates the collapse of larger 
segments of molecular hygrogen clouds than would 
occur without this interaction, resulting in higher- 
mass stellar associations. 

8. Even with the resolution of Hubble ACS/WFC, 
main sequence turnoffs are not evident for the vast 
majority of stellar associations in M51, due to the 
limiting magnitude of our images, and the vertical 
nature of the turnoff at the bright end of the main 
sequence. 

This type of study has not been possible before because 
no other face-on spiral has been completely imaged in 
the photometric bands necessary to sufficient depth and 
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resolution to conduct such an investigation on an entire 
external galaxy. In recognition of this fact, the study 
itself has been very challenging and while we have had 
some success in measuring the ages of associations in 
our sample, there have been many fundamental limits to 
what we can achieve and access because of the distance 
to the galaxy, the limiting magnitude of the datasets, 
and the available band widths, even when done with the 
exquisite optics of the HST from space. 

Despite these limitations, we do believe the dataset 
provides an unprecedented view of the very recent SF 



process across the entire disk of the galaxy and reflects a 
definitive change in the SFR from one side of the galaxy 
to the other probably as a direct result of the interaction 
with NGC 5195. This paper represents only the first 
part of a larger study to study the various SF complexes 
in M51 using multiwavelength datasets including addi- 
tional data from GALEX and Spitzer. That work will be 
reported in subsequent publications. 

We acknowledge financial support of this work by 
NASA/STScI Grant HST-AR-10684.01-A. 
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TABLE 1 

Stellar Association Parameters for 120 Single- Aged Associations in the M51/NGC 5195 System 



Identifier 


RA a 




Dcc. a 


M Bi»t C 


M v int c 


M /i»t C 


Ms/ 




M/, d 


Radius 
(pc) e 


Mass 
(M ) 


Age 
(Myrs) 


Cl022x3960yl058x4062yll47 


13:30:6.6247 


47: 


14:33.358 


-9.42 


-9.51 


-10.0 


-8.19 


-8.25 


-8.55 


61.0 


1.37xl0 5 


86 


c3021x450y731x488y760 


13:29:58.297 


47: 


14:11.042 


-10.2 


-10.3 


-11.0 


-9.64 


-9.82 


-10.4 


29.4 


5.05xl0 5 


200 


c3021x489y601x526y630 


13:29:58.897 


47: 


14:9.6002 


-10.0 


-10.1 


-10.5 


-8.91 


-9.05 


-9.23 


21.4 


2.41x10 s 


86 


c3021x708yl939x758yl988 


13:29:52.393 


47: 


13:55.922 


-10.0 


-10.1 


-10.5 


-8.97 


-9.11 


-9.47 


36.0 


2.47x10 s 


86 


c3021x782yl340x816yl376 


13:29:55.294 


47: 


13:53.766 


-8.93 


-9.00 


-9.41 


-7.93 


-8.03 


-8.19 


26.4 


8.43xl0 4 


86 


c3021x889y657x923y688 


13:29:58.535 


47: 


13:49.797 


-10.6 


-10.5 


-10.4 


-9.89 


-9.81 


-9.69 


24.1 


2.34xl0 4 


4 


c3021xll24yl635xll43yl662 


13:29:53.785 


47: 


13:36.833 


-9.12 


-9.00 


-8.96 


-8.24 


-8.22 


-8.30 


22.8 


5.61x10 s 


4 


c3021xll54yl610xll77yl651 


13:29:53.902 


47: 


13:35.405 


-8.24 


-8.22 


-8.76 


-7.17 


-7.24 


-7.88 


33.5 


3.31 xlO 4 


57 


c3021xl508y293xl566y332 


13:30:0.0732 


47: 


13:19.557 


-9.72 


-9.80 


-10.3 


-8.69 


-8.94 


-9.56 


39.7 


2.40xl0 5 


140 


c3021xl533yl36xl596yl90 


13:30:0.7910 


47: 


13:18.115 


-10.5 


-10.4 


-10.9 


-9.95 


-9.97 


-10.2 


55.0 


2.59xl0 4 


5 


c3021x2008yl044x2042yl071 


13:29:56.400 


47: 


12:54.000 


-8.28 


-8.21 


-8.19 


-7.66 


-7.59 


-7.38 


25.0 


2.62xl0 3 


4 


c3021x2016yl329x2042yl352 


13:29:55.056 


47: 


12:53.286 


-8.22 


-8.21 


-8.57 


-7.59 


-7.61 


-7.71 


13.8 


3.09xl0 3 


5 


c3021x2031yl397x2057yl430 


13:29:54.697 


47: 


12:52.558 


-9.07 


-9.05 


-9.35 


-8.11 


-8.13 


-8.12 


18.5 


6.68xl0 3 


5 


c3021x2122ylll8x2145yll41 


13:29:56.015 


47: 


12:48.603 


-8.06 


-8.00 


-8.09 


-7.62 


-7.54 


-7.47 


14.2 


2.17x10 s 


4 


c3021x2847yl871x2878yl907 


13:29:52.199 


47: 


12:11.524 


-9.53 


-9.52 


-9.48 


-8.82 


-8.82 


-8.68 


23.0 


8.36x10 s 


4 


c3021x3532y20x3577y56 


13:30:0.8642 


47: 


11:40.556 


-10.4 


-10.4 


-10.7 


-9.38 


-9.37 


-9.56 


30.3 


2.40xl0 4 


5 


c3021x3823yl487x3857yl521 


13:29:53.833 


47: 


11:23.995 


-8.30 


-8.19 


-8.65 


-7.19 


-7.31 


-7.65 


27.7 


3.29x10 s 


5 


c3021x3842yl462x3880yl493 


13:29:53.928 


47: 


11:22.923 


-8.15 


-8.21 


-8.92 


-6.67 


-6.81 


-7.04 


36.7 


4.91 xlO 4 


110 


c3021x3898y716x3935y757 


13:29:57.454 


47: 


11:21.481 


-8.95 


-9.00 


-9.37 


-8.10 


-8.18 


-8.29 


31.0 


6.20x10 s 


5 


c3022x517yl929x543yl960 


13:30:2.5195 


47: 


14:10.685 


-7.26 


-7.32 


-7.77 


-5.82 


-6.06 


-6.38 


14.5 


1.82xl0 4 


86 


c3022x628yl633x683yl680 


13:30:3.9111 


47: 


14:4.9172 


-7.89 


-7.85 


-8.13 


-6.51 


-6.53 


-6.54 


24.7 


2.24x10 s 


5 


c3022x904yl815x949yl860 


13:30:2.9772 


47: 


13:50.525 


-8.04 


-8.15 


-8.74 


-6.49 


-6.72 


-6.93 


28.1 


5.59xl0 4 


160 


c3022xl022yl8291058yl866 


13:30:2.8784 


47: 


13:45.114 


-8.55 


-8.61 


-9.15 


-7.58 


-7.64 


-7.98 


22.5 


6.08xl0 4 


86 


c3022xll86y571xl205y590 


13:30:9.1186 


47: 


13:39.360 


-7.16 


-7.14 


-7.41 


-6.72 


-6.76 


-6.73 


20.3 


1.14x10 s 


5 


c3022xl435yl792xl470yl826 


13:30:2.9992 


47: 


13:24.597 


-7.93 


-8.09 


-8.59 


-6.88 


-7.00 


-7.25 


22.6 


5.05xl0 4 


160 


c3022xl521yl530xl555yl564 


13:30:4.2736 


47: 


13:20.999 


-8.86 


-8.75 


-8.86 


-8.05 


-8.05 


-8.06 


31.9 


4.49x10 s 


4 


c3022xl633y721xl664y746 


13:30:8.2800 


47: 


:13:16.673 


-7.87 


-8.18 


-8.86 


-7.25 


-7.52 


-7.80 


32.3 


8.80X10 4 


310 


c3022xl680yl247xl705yl277 


13:30:5.6396 


47: 


13:13.446 


-8.83 


-8.85 


-9.43 


-8.04 


-8.15 


-8.71 


15.3 


7.82 xlO 4 


86 


c3022xl708y820xl723y844 


13:30:7.7526 


47: 


13:13.075 


-8.28 


-8.29 


-9.20 


-7.88 


-7.85 


-8.42 


16.2 


2.15xl0 4 


22 


c3022xl842yl441xl868yl463 


13:30:4.6801 


47: 


13:5.1654 


-9.26 


-9.21 


-9.77 


-8.23 


-8.20 


-8.57 


11.1 


8.27x10 s 


5 


c3022xl875yl444xl915yl469 


13:30:4.6545 


47: 


13:2.9956 


-8.78 


-8.87 


-9.43 


-8.01 


-8.14 


-8.31 


26.7 


7.66X10 4 


86 


c3022xl965yl585x2005yl614 


13:30:3.9367 


47: 


12:58.683 


-8.08 


-8.18 


-9.05 


-6.92 


-7.10 


-7.79 


17.0 


8.30X10 4 


250 


c3022x2085yll65x2106yll88 


13:30:5.9985 


47: 


12:53.643 


-7.90 


-7.93 


-8.37 


-7.42 


-7.42 


-7.65 


28.0 


2.37x10 s 


5 


c3022x2101y588x2126y614 


13:30:8.8073 


47: 


12:53.286 


-7.50 


-7.57 


-7.96 


-6.75 


-6.80 


-6.95 


16.9 


1.65x10 s 


5 


c3022x2309yl481x2332yl505 


13:30:4.3688 


47: 


12:42.121 


-7.13 


-7.12 


-7.30 


-6.45 


-6.53 


-6.72 


27.4 


1.11x10 s 


5 


c3022x2324yl687x2354yl731 


13:30:3.3105 


47: 


12:40.679 


-8.31 


-8.37 


-8.94 


-7.38 


-7.55 


-8.00 


27.2 


4.91 xlO 4 


86 


c3022x2449yl572x2491yl611 


13:30:3.8891 


47: 


12:34.554 


-9.71 


-9.85 


-10.2 


-9.01 


-9.22 


-9.47 


43.7 


1.76x10 s 


86 


c3022x2480yl085x2514ylll2 


13:30:6.2878 


47: 


12:33.483 


-7.65 


-7.90 


-8.83 


-6.54 


-6.84 


-7.65 


32.1 


9.92 xlO 4 


430 


c3022x2618yl453x2655yl485 


13:30:4.4384 


47: 


12:26.273 


-10.0 


-9.99 


-10.0 


-9.41 


-9.37 


-9.36 


37.1 


1.39X10 4 


4 


c3022x2629y959x2659y986 


13:30:6.8884 


47: 


12:26.644 


-9.72 


-9.63 


-9.63 


-8.99 


-8.89 


-8.71 


21.8 


9.92x10 s 


4 


c3022x2709yl546x2723yl568 


13:30:3.9843 


47: 


12:22.318 


-7.50 


-7.41 


-7.75 


-7.12 


-7.03 


-7.31 


5.40 


1.54x10 s 


5 


c3022x2735yl709x2753yl728 


13:30:3.1677 


47: 


12:20.519 


-7.27 


-7.33 


-7.76 


-6.60 


-6.68 


-7.07 


17.6 


1.84 xlO 4 


86 


c3022x3101yl813x3143yl850 


13:30:2.5671 


47: 


12:1.4447 


-8.82 


-8.92 


-9.35 


-8.08 


-8.22 


-8.52 


40.0 


7.75 xlO 4 


86 


c3022x3561y786x3605y846 


13:30:7.4890 


47: 


11:39.842 


-8.96 


-9.28 


-10.3 


-8.08 


-8.51 


-9.55 


47.5 


4.34 x10 s 


540 


c3022x3861yl479x3894yl527 


13:30:3.9843 


47: 


11:24.365 


-7.62 


-7.78 


-8.72 


-6.58 


-6.67 


-7.51 


25.2 


8.43X10 4 


390 


c3022x3861y936x3908y973 


13:30:6.7456 


47: 


11:24.722 


-9.51 


-9.47 


-9.43 


-8.79 


-8.77 


-8.68 


33.0 


8.26x10 s 


4 


c402 lxl 738y32xl 775y74 


13:29:41.953 


47: 


13:1.1965 


-8.01 


-8.04 


-8.70 


-7.20 


-7.30 


-7.38 


33.0 


3.96X10 4 


96 


c4021x2677y615x2717y650 


13:29:38.902 


47: 


12:13.680 


-8.69 


-8.72 


-8.91 


-8.14 


-8.17 


-8.19 


25.5 


4.70x10 s 


5 


c4021x2711y385x2759y443 


13:29:39.960 


47: 


12:12.238 


-9.14 


-9.27 


-9.67 


-8.60 


-8.73 


-8.99 


41.8 


1.04 x10 s 


86 


c4021x3433y239x3487y286 


13:29:40.462 


47: 


11:36.958 


-8.08 


-8.28 


-9.43 


-7.36 


-7.63 


-8.64 


54.4 


6.20x10 s 


7 


c4021x3464y32x3497y67 


13:29:41.495 


47: 


11:36.244 


-9.13 


-9.19 


-9.81 


-8.34 


-8.44 


-8.78 


20.0 


9.77X10 4 


77 


c4021x3501y776x3544y818 


13:29:37.921 


47: 


11:33.003 


-8.80 


-8.90 


-9.59 


-7.45 


-7.64 


-8.15 


28.6 


9.79X10 4 


130 


c4021x3707y534x3744y571 


13:29:39.023 


47: 


11:23.280 


-8.00 


-7.85 


-7.91 


-7.16 


-7.06 


-6.90 


22.8 


2.00x 10 s 


4 


c4022xl86yl69x217yl97 


13:29:51.936 


47: 


14:21.479 


-7.90 


-7.94 


-8.30 


-7.17 


-7.19 


-7.16 


25.5 


2.34x10 s 


5 


c4022x587y467x642y497 


13:29:50.375 


47: 


14:0.6051 


-9.06 


-9.13 


-9.28 


-8.26 


-8.31 


-8.23 


44.4 


6.65x10 s 


5 


c4022x692yl002x744yl066 


13:29:47.640 


47: 


13:54.480 


-8.65 


-8.75 


-9.25 


-7.54 


-7.71 


-7.93 


55.6 


6.72X10 4 


86 


c4022x873y435x910y474 


13:29:50.471 


47: 


13:46.556 


-8.00 


-8.50 


-9.42 


-6.90 


-7.45 


-8.26 


36.3 


2.07x10 s 


610 


c4022x880y668x911y698 


13:29:49.321 


47: 


13:46.199 


-8.20 


-8.27 


-8.61 


-7.59 


-7.77 


-7.98 


19.9 


3.10x10 s 


5 


c4022xl032yl210xl081yl238 


13:29:46.633 


47: 


13:37.561 


-9.44 


-9.39 


-9.48 


-8.55 


-8.52 


-8.57 


33.9 


7.79x10 s 


4 


c4022xl093yll86xlll5yl208 


13:29:46.750 


47: 


13:35.034 


-8.94 


-9.07 


-9.58 


-8.21 


-8.48 


-8.80 


22.4 


1.27x10 s 


160 


c4022xll58y346xll80y376 


13:29:50.855 


47: 


13:32.878 


-7.71 


-7.88 


-8.79 


-6.73 


-6.88 


-7.44 


25.7 


9.14X10 4 


390 


c4022xl713yll8xl738yl61 


13:29:51.793 


47: 


13:5.5224 


-8.72 


-8.66 


-8.76 


-8.41 


-8.37 


-8.32 


22.6 


3.98x10 s 


4 


c4022xl775yl05xl807yl39 


13:29:51.936 


47: 


13:2.2814 


-8.08 


-8.15 


-9.04 


-7.32 


-7.44 


-7.68 


35.4 


8.22X10 4 


250 
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Identifier 




RA a 




Dcc. a 


M Bi»t C 


M v int c 


M /i»t C 


Ms/ 


M 


v, d 


M/, d 


Radius 




Mass 


Age 




























(pc) e 


(M ) 


(Myrs) 


c4022xl965y456xl995y487 


13 


29: 


:50.185 


47: 


12:52.558 


-8.48 


-8.44 


-8.39 


-8.00 


-7. 


96 


-7.78 


zo . 


,o 


Q 
O 


1 7 v 1 n3 

. 1 t X 1U 


4 


c4022x2 109y 746x2 1 56y 785 


13: 


29: 


:48.647 


47: 


12:44.277 


-8.09 


-8.23 


-9.03 


-7.41 


-7 


.63 


-8.02 


A 1 


.0 


7 


.0 i X 1U 


220 


c4022x2149yl043x2199yl084 


13: 


:29: 


47.182 


47: 


12:41.764 


-8.56 


-8.76 


-9.54 


-7.94 


-8. 


.16 


-8.62 


oo. 


A 
. i 


1 


.OZ X _LU 


250 


c4022x2289yl064x2329yll03 


13: 


:29: 


:47.065 


47: 


12:35.282 


-8.01 


-8.15 


-8.87 


-7.01 


-7. 


.21 


-7.63 


ol. 


Q 

,o 


D 


-4oX _LU 


200 


c4022x2457yl512x2491yl553 


13: 


'29: 


: 44.831 


47: 


12:26.273 


-9.61 


-9.58 


-9.74 


-8.81 


-8 


.77 


-8.71 


oo. 


,Z 


O 

y 


,1 ( X IU 


4 


c4022x2554y502x2612y566 


13: 


29: 


:49.727 


47: 


12:22.318 


-9.34 


-9.39 


-10.0 


-8.59 


-8. 


.63 


-9.00 


36. 


,4 


i 
i 


,40 X ±U 


110 


c4022x2787y25x2821y58 


13: 


:29: 


52.104 


47: 


12:12.238 


-9.43 


-9.40 


-9.35 


-8.78 


-8. 


76 


-8.67 


24. 


5 


( 


f^o \/ 1 n3 
,oy x iu 


4 


c4022x2826yl550x2879yl603 


13 


29: 


: 44.520 


47: 


12:7.1987 


-7.99 


-8.03 


-8.84 


-6.64 


-6 


70 


-7.03 


50. 


,7 


1 


,yo x iu 


31 


c4022x2840y571x2883y602 


13: 


'29: 


:49.390 


47: 


12:8.2836 


-8.94 


-9.03 


-9.63 


-8.33 


-8. 


55 


-9.15 


32 


,7 


1 


,loX IU 


140 


c4022x2945y717x2974y749 


13 


29: 


: 48.673 


47: 


12:3.2437 


-8.12 


-8.04 


-8.74 


-7.48 


-7. 


49 


-7.79 


13. 


,8 


z 


,oy x iu 


40 


c4022x3333yl755x3350yl779 


13 


'29: 


: 43.487 


47: 


11:43.083 


-8.26 


-8.29 


-8.57 


-7.40 


-7. 


50 


-7.76 


11. 


,7 





-Zl X IU 


5 


c4022x3437yl043x3468yl084 


13 


29: 


:46.944 


47: 


11:38.400 


-7.96 


-7.82 


-8.59 


-7.48 


-7. 


37 


-7.53 


38. 


,6 


1 


.O i X IU 


24 


c4022x3868y340x3912y396 


13: 


'29: 


:50.302 


47: 


11:17.526 


-9.25 


-9.36 


-10.0 


-8.44 


-8. 


.67 


-9.30 


46. 


.2 


1 

1 


.00 X IU 


180 


c502 Ix264y445x3 14y509 


13: 


'29: 


58.967 


47: 


11:5.6341 


-11.3 


-11.2 


-11.4 


-10.5 


-10.5 


-10.6 


41. 


,6 


A 
1 


t^fi \y 1 n4 


4 


c5021x665yl627x700yl662 


13: 


29: 


:53.280 


47: 


10:44.046 


-10.8 


-10.9 


-11.1 


-10.0 


-10.1 


-10.3 


24. 


3 


O 


,4y x ±u 


5 


c5021x669yl626x701yl667 


13: 


:29: 


53.280 


47: 


10:44.046 


-10.8 


-10.9 


-11.1 


-10.2 


-10.2 


-10.4 


32. 


,0 


O 


,01 X ±u 


5 


c5021x703y849x726y885 


13: 


:29: 


:56.975 


47: 


10:43.675 


-7.73 


-8.00 


-9.13 


-6.63 


-7 


11 


-8.02 


21. 


,0 


1 
1 


-4U X1U 


540 


c5021x739y 1 142x770yl 180 


13: 


'29: 


:55.583 


47: 


10:41.519 


-8.33 


-8.32 


-8.86 


-7.76 


-7. 


.80 


-7.94 


29. 


1 


O 


ci v 1 n4 
,Dl X IU 


57 


c5021x791yl253x813yl271 


13: 


29: 


:55.081 


47: 


10:38.635 


-7.31 


-7.44 


-8.16 


-6.50 


-6 


67 


-6.98 


9.30 


Q 
O 


.00X 11) 


200 


c5021x791yl620x828yl644 


13 


'29: 


53.327 


47: 


10:37.921 


-8.88 


-8.96 


-9.49 


-7.93 


-8. 


19 


-8.66 


34 


5 


O 
O 


07vi ri4 

.Z 1 X IU 


86 


c5021xll41yl883xll67yl934 


13 


29: 


:51.936 


47: 


10:20.274 


-8.49 


-8.49 


-9.10 


-7.40 


-7. 


43 


-7.80 


37 


2 


A 
y l 


oc v 1 fl4 
,Z0 X IU 


57 


c5021xll49yl574xll82yl604 


13: 


'29: 


:53.422 


47: 


10:20.274 


-9.58 


-9.58 


-10.1 


-8.72 


-8. 


75 


-9.06 


29. 


.0 


1 
1 




57 


c5021xll69yl324xl217yl377 


13 


29: 


:54.576 


47: 


10:19.560 


-10.5 


-10.5 


-11.1 


-9.68 


-9 


77 


-10.2 


35. 


,8 


9 
Z 


.OO X IU 


57 


c5021xll78yl641xl213yl673 


13: 


:29: 


53.111 


47: 


10:18.846 


-9.10 


-9.07 


-9.77 


-8.32 


-8. 


.34 


-8.78 


37. 


,8 





70 sy 1 ri4 


38 


c5021xl253yl509xl284yl540 


13: 


:29: 


:53.712 


47: 


10:15.605 


-9.06 


-8.90 


-8.73 


-8.44 


-8. 


.33 


-8.03 


28. 


,8 





on \y 1 n3 
-ZUX IU 


4 


c5021xl735yl344xl775yl378 


13: 


'29: 


:54.360 


47: 


9:51.833 


-8.35 


-8.41 


-8.87 


-7.42 


-7 


.47 


-7.66 


35. 


,6 


A 
1 


07 \/ 1 n4 
.y i X IU 


86 


c5021xl852y874xl899y907 


13: 


'29: 


56.568 


47: 


9:46.793 


-9.22 


-9.33 


-9.86 


-8.60 


-8. 


.67 


-8.98 


32. 


3 


1 

1 


-DoX IU 


160 


c5021x2324yl077x2369yll34 


13: 


29: 


:55.440 


47: 


9:23.035 


-8.62 


-8.69 


-9.32 


-7.75 


-7. 


.78 


-8.01 


50. 





7 


nn \/ 1 n4 


96 


c5021x2429yl611x2449yl636 


13: 


'29: 


52.943 


47: 


9:17.638 


-8.15 


-8.02 


-7.95 


-7.66 


-7. 


57 


-7.36 


14. 


3 


Z 


.ZoX iU 


4 


c5021x2430yl21x2470yl57 


13: 


30: 


0.0256 


47: 


.9:19.794 


-8.83 


-8.83 


-9.07 


-8.04 


-8. 


10 


-8.19 


35. 


3 





qqv 1 n3 
.00X IU 


5 


c5021x2537y414x2560y438 


13: 


'29: 


:58.630 


47: 


9:14.397 


-8.61 


-8.60 


-8.79 


-8.02 


-8. 


05 


-8.23 


22. 


6 


A 
y l 


o7 V 1 n3 


5 


c5021x2552y318x2590y354 


13 


29: 


59.040 


47: 


9:13.326 


-8.90 


-8.90 


-9.14 


-8.03 


-8. 


04 


-8.11 


34. 


,4 





kf. \/ 1 n3 
,DDX IU 


5 


c5021x2806yl020x2840yl051 


13: 


:29: 


55.631 


47: 


9:0.0054 


-8.51 


-8.36 


-8.31 


-7.80 


-7. 


.67 


-7.41 


26. 


,5 


O 


1 S v 1 f|3 


4 


c5021x2854y667x2890y702 


13 


29: 


:57.286 


47: 


8:57.835 


-9.73 


-9.66 


-9.61 


-9.28 


-9. 


,24 


-9.13 


25. 





1 


nn \/ 1 n4 

.UUX IU 


4 


c5021x3033y439x3063y481 


13: 


'29: 


:58.344 


47: 


8:49.925 


-8.38 


-8.39 


-8.64 


-7.72 


-7. 


.78 


-7.99 


43. 


,4 


O 
O 


-OOX IU 


5 


c5021x3054yl592x3074yl629 


13: 


:29: 


52.822 


47: 


8:47.041 


-8.54 
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a Epoch J2000 

k All magnitudes in the Vega magnitude system 

c Total integrated magnitude from a circular aperture of radius and central position listed in Table 1. 

^ Total magnitude from summing all sources in the association. This value will not include diffuse light from unresolved stars, and is thus a lower 
limit of the magnitude of the association. 

e The average radius of the association as determined by the average FWHM of the flux density distribution. 



